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Abstract 

We study the effect of weak localization on the transition temperatures of 
superconductors using time-reversed scattered state pairs, and find that the 
weak localization effect weakens electron-phonon interactions. With solving 
the BCS T c equation, the calculated values for T c are in good agreement with 
experimental data for various two- and three-dimensional disordered super- 
conductors. We also find that the critical sheet resistance for the suppression 
of superconductivity in thin films does not satisfy the universal behavior but 
depends on sample, in good agreement with experiments. 



PACS numbers: 74.20.-z, 74.40.+k, 74.60.Mj 
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I. INTRODUCTION 



Since the scaling theory of Anderson localization, 1 our understanding of the electronic 
properties of disordered conductors has been advanced considerably. 2 ' 3 However, the effect 
of localization on superconductivity is still not well understood. 4-7 A unified picture for 
the disorder effects on the superconducting temperature (T c ) and the critical field (H c2 ) is 
still lacking. In the presence of disorder, it was even claimed that the Eliashberg theory 
breaks down for two-dimensional superconductors. 8 There have been many experimental 
studies to explain a competition between localization and superconductivity in two- and 
three-dimensional systems. 8-16 In homogeneous amorphous thin films, an empirical formula 5 
showed that the reduction of T c is proportional to the sheet resistance R n - For bulk amor- 
phous InO x , Fiory and Hebard 13 found that both the normal-state conductivity a and the 
critical temperature vary as (k F £)~ 2 due to the localization effect, where kp and t are the 
Fermi wave vector and the elastic mean free path, respectively. In A15 superconductors 
such as Nb 3 Sn, V 3 Si, and V 3 Ga, the universal degradation of T c with impurities was also 
found, 17 following the variation of p 2 oc (k F £)~ 2 , 18-20 where p is the normal-state resistivity. 
Previously, the decrease of T c with increasing of disorder was attributed to the enhanced 
Coulomb repulsion yU*. 4 ' 17 ' 21 ' 22 However, tunneling experiments 8 ' 23-25 did not support such 
an argument and indicated instead a decrease of the electron-phonon coupling A. 8,23-26 As 
an alternative explanation, a singularity of the density of states at the Fermi level caused by 
long-range Coulomb interactions was suggested to change T c 's for Pb and Sn thin films. 27 
However, it was pointed out that the singularity of the density of states does not affect 
thermodynamic quantities such as T c . 28 In addition, it was shown experimentally that the 
effect of the Coulomb gap is minor in superconductors. 29 ' 30 

In this paper we present the results of theoretical studies on the weak localization effect 
on the superconducting temperatures of disordered superconductors. Using time-reversed 
scattered state pairs, we are able to explain available experimental data, based on a uni- 
fied picture that the electron-phonon coupling decreases with increasing of disorder. Our 
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theory also explains other experimental results that in disordered superconductors T c can 
be enhanced by spin-orbit scatterings, 12,31-33 which results from the anti-localization effect. 
In thin films, we suggest that the critical sheet resistance for the suppression of supercon- 
ductivity is not a universal constant but a sample-dependent quantity. Some preliminary 
results were reported elsewhere. 34 



II. THEORY 

It has been recently realized that there exist localization corrections in electron-phonon 
interactions. To find these correction terms, 35-38 it is a prerequisite to understand the limita- 
tion of Anderson's theorem and the pairing problem in Gor'kov's formalism and Bogoliubov- 
de Gennes equations. The Anderson's theory in dirty superconductors 39 is based on the fact 
that the exact eigenstates in the presence of impurities consist of time-reversed degenerate 
pairs; the scattered state ip n of an electron with spin up is paired with the other electron with 
spin down in the time-reversed state ip n - Then, the reduced Hamiltonian in scattered-state 
representation is written as 

H 'red = 2 V nn'clcl,C n C n , (1) 
nn' 

where and c n are creation and annihilation operators, respectively, for an electron in the 
state ip n . Assuming a point coupling —V5{r 1 — r 2 ) for phonon-mediated electron-electron 
interactions, the matrix elements V nn > are expressed as 40 

V nn , = -V j <,(r)^,(rty fi (rty n (r)dr. (2) 

If the scattered state ip na is expanded in terms of plane waves 

Vw = Y.h*$\n), (3) 

k 

V nn r is rewritten as 41 

V m , = -V(l+ £ (-k'\n)(k\nr(k-q\n')(k-q\n')(-k' -?>')*)■ ( 4 ) 
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The second term in Eq. (4) is negligible in dirty limit, where the mean free path I is in the 
range of 100 A, while its contribution is meaningful in weak localization limit, where I is of 
the order of 10 A. Consequently, the Anderson's theorem is only valid in the dirty limit. 

In order to calculate the correction term in Eq. (4), we need information on the weakly 
localized scattered states. Kaveh and Mott 42 derived these scattered states in the form 
of power-law and extended wavefunctions for both two and three dimensional systems. 
Haydock 43 also showed that the asymptotic form of the scattered states is power-law like 
in weakly disordered two-dimensional systems. Here we use the scattered states of Kaveh 
and Mott 42 to calculate the matrix elements V nn i. In this case, since we are dealing with 
the bound state of a Cooper pair in a BCS condensate, only the power-law wavefunctions 
within the BCS coherence length are relevant. 44 In weak localization limit, the effective 
coherence length is reduced to £ e // ~ \fWo- It is important to take into account the size of 
the Cooper pairs for self-consistent calculations of the matrix elements. A similar situation 
also occurs for the localized states, which are not very sensitive to the change of the bound- 
ary conditions. 45 Because of the power-law-like wavefunctions, the correction term due to 
weak localization represents physically the decrease of the amplitudes of the plane waves. 
Thus, the Cooper-pair wavefunctions basically consist of the plane waves with reduced am- 
plitudes. Including the weak localization effect, the resulting matrix elements V nn > for two- 
and three-dimensional systems are written as 

V% = -V[l - -2-ln(L/e)], (5) 

Because of the impurity effect, the electron-phonon interactions are weakened, as clearly 
seen in Eqs. (5) and (6). 
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III. RESULTS AND DISCUSSION 



If we assume that the electron-phonon coupling constant is modified by the weak lo- 
calization effect, the effective coupling constant X e ff from Eq. (5) for three-dimensional 
systems can be written as 



Kjj = A[l - 




Then, using the modified BCS T c equation, 

k B T c = 1.13hu D exp(-l/\ e ff), (8) 

the change of T c relative to T co (for a pure metal) can be easily estimated to first order in 
the weak localization correction term of the coupling constant; 

For a metal with oo D = 300 K, k F = 1 A -1 , and L = 1000 A/T, the variations of T c 
are plotted as a function of l/k F £ in Fig. 1, assuming T co = 4 and 12 K. The inelastic 
mean free path results from electron-electron scatterings, and the correction term in Eq. 
(9) is negligible for I»f. We find that T c changes slowly with increasing of l/k F £ until 
\jk F i equals to 0.1, in good agreement with experimental results. 13 ' 18-20 This behavior which 
satisfies the Anderson's theorem is attributed to the fact that the change of T c is proportional 
to p 2 , i.e., the square of impurity concentration. However, for l/k F £ > 0.1, T c decreases 
more significantly due to the weak localization effect caused by ordinary impurities. In this 
weak localization limit, the Anderson's theorem is not valid. Our theoretical results are also 
compared with experimental data 46 for A15 superconducting materials in Fig. 2. Using the 
same values of k F = 0.87 A -1 and L = 1000 A/T for Nb 3 Ge and V 3 Si with u D = 302 and 
330 K and T co = 23 and 17 K, respectively, we find good agreements between theory and 
experiment. In this case, since it is difficult to evaluate k F £ up to a factor of 2, 47 we assume 
that p = lOOpQcm corresponds to k F £ = 3.3 and 3.75 for Nb 3 Ge and V 3 Si, respectively. 
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Similarly, we can write down the effective electron-phonon coupling constant for two- 
dimensional systems such as 

\ eff = \[l--^- l ln{L/i)]. (10) 
Since l/k F £ is related to R n by the Drude formula, the change of T c is expressed as, 



T C o T c 1 6 



2 



Xn 2 h 



R n ln(L/£), (11) 



and this formula well satisfies the empirical relationship between T c and R n for two- 
dimensional superconductors. In Fig. 3, the variations of T c with l/kp£ are drawn for 
two superconductors with T co = 4 and 8 K, assuming the same values of uo D = 300 K, i 
= 4 A, and L = 1000 ky/T. The inelastic mean free path obtained from disordered two- 



dimensional systems 48 are employed, with the 1/ \J InT dependence removed. In contrast 
to the 3-dimensional case, T c varies linearly with increasing of impurity concentration, and 
the initial slope of T c depends on superconductor because of the prefactor 1/A in Eq. (12). 
Thus, the critical sheet resistances for the suppression of superconductivity in thin films do 
not provide the universal behavior, in good agreement with experiments. 5 ' 9,12 

Since the inelastic mean free path depends on temperature, the effective electron-phonon 
coupling also varies with temperature. Including the temperature effect on X e ff, the variation 
of the gap parameter with temperature is plotted for two different values of l/k F £ = 0.039 
and 0.078 in Fig. 4. The gap parameters are found to be lower than those obtained using 
the T-independent A e //'s. Experimentally, this behavior may be observable for lower T c 
superconductors, while strongly-coupled superconductors such as Pb and Nb may not be 
appropriate because the BCS model is not applicable for these systems. 

For Mo-C 49 and a-MoGe 12 thin films, our calculated T c 's are plotted as a function of 
l/kp£ and compared with experimental data in Fig. 5. Here we employ the Drude formula 
to represent the decrease of £ when the sheet resistance R n increases. As in the 3-dimensional 
case, because of the difficulty in evaluating l/k F £ up to a factor of 2 from experimental data, 
we assume I = 3.5 A and L = 2000 kVT for best fit in the Mo-C sample. For the a-MoGe 
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film, the use of i = 4.0 Aand L = 1000 ky/T is found to give the best agreement with 
experiment. In this case, the Drude formula is slightly adjusted by the relation = 
1.6667 (e 2 /2nh)R n . Similarly to the three-dimensional case, we find that the critical sheet 
resistance for the suppression of superconductivity does not follow the universal behavior. 

We can also examine the effect of weak localization on superconductivity, using the strong 
or weak coupling Green's function theory. In previous approaches, 50 " 52 because the Ander- 
son's time-reversed scattered-state pairs were not employed, the correction term due to weak 
localization was missing in the electron-phonon interaction. The Green's function formalism 
leads to the pairing states formed by a linear combination of the scattered states. 35-38 Since 
these are still extended states, we do not expect the weak localization effect. Using the 
Anderson's pairing condition for the strong coupling equation and the Einstein model for 
the phonon spectrum, the gap equation can be written as, 37 

A(n, u) = £ X(co - co') £ V nn ,^^±, (12) 

where 

V„ n > = -vJ |i/.„(r)| 2 |i*„.(r)| 2 rfr, (13) 

In this case, the use of V nn > in Eq. (2) can also give rise to a modification of the electron- 
phonon interaction due to impurities in the strong coupling theory. In previous theories, 50-52 
however, the electron-phonon interaction remains unchanged, even if the wavefunctions are 
localized. 

It is interesting to see the superconductor-to-insulator transitions in ultrathin films, 
which usually occur by changing film thickness or applying magnetic fields. 49,53 ' 54 Previously, 
a dirty boson theory 55 was used to explain this experimental feature, implying that the 
critical sheet resistance for the suppression of supercondcutivity is a universal constant, 
h/4e 2 = 6.45 KQ. However, although it is difficult to determine experimentally a well- 
defined value for the critical sheet resistance because the transitions do not occur sharply, 



the measured values for were shown to depend on sample. 49 ' 53 ' 54 We point out that the 
weak localization effect considered here is still a dominant contribution to the suppression 
of superconductivity over other higher order terms beyond the weak localization regime. 2,3 
In fact, Fig. 3 shows that T c is completely suppressed in the region of kpi ss 6 or 7, which 
is still in the weak localization regime. Thus, our results indicate that the superconductor- 
insulator transition is not a sharp phase transition but a crossover phenomena from quasi-two 
dimensional to two-dimensional. 

Finally, we suggest that if the decrease of T c is caused by weak localization in disordered 
superconductors, adding impurities with large spin-orbit couplings will compensate for this 
decrease. In fact, such behavior was observed for several 3-dimensional samples, 12 ' 31 ~ 33 while 
it needs to be tested for 2-dimensional case. We expect that the critical sheet resistance is 
increased by enhancing the spin-orbit scattering. It is known that magnetic fields suppresses 
the weak localization effect. In this case, however, since the magnetic field decreases T c , the 
T c decrease in a pure sample should be compared with that of a weakly disordered one. If 
the difference would exist, it is suggested to result from the weak localization effect, then, 
the critical sheet resistance will also change with increasing of magnetic field. 

IV. CONCLUSIONS 

In conclusion, we have studied the effect of weak localization on superconductors within 
the BCS theory. We find that the weak localization decreases the electron-phonon coupling 
constant, thereby, suppressing T c . The calculated variations of T c with increasing of im- 
purity concentration are found to be in good agreement with experiments for both 2- and 
3-dimensional systems. The recovery of T c with impurities having large spin-orbit scatterings 
supports strongly our theory. We suggest that the critical sheet resistance for the suppres- 
sion of superconductivity in thin films is not a universal constant, but a sample-dependent 
quantity, in good agreement with experiments. 
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FIGURES 

FIG. 1. Variation of T c with disorder parameter \jkpt (which represents ordinary impurity 
concentration) for 3-dimensional superconductors with T co = 4 and 12 K. 

FIG. 2. Calculated T c 's versus resistivity p for 3-dimensional Nh>3Ge (dotted line) and VsSi 
(solid line). Experimental data are from Ref. 46. 

FIG. 3. Variation of T c with disorder parameter l/fcp£ for 2-dimensional superconductors with 
T co = 4 and 8 K. 

FIG. 4. Temperature dependence of the gap parameter for 2-dimensional superconductors. 
The values of \jkpt = 0.039 and Rn = 2000 O are chosen for the upper curves, while l/fc^ = 0.078 
and Rn = 3000 Q for the lower curves, with T co = 4 K fixed. The T-independent (T-dependent) 
effective electron-phonon coupling is used for the solid (dotted) lines. 

FIG. 5. Calculated T c 's versus sheet resistance Rn for a-MoGe (solid line) and Mo-C (dotted 
line) thin films. Experimental data for a-MoGe and Mo-C are from Refs. 12 and 49, respectively. 
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